Exciton-Polariton scattering for defect detection in cold atom Optical Lattices 



Hashem Zoubi, and Helmut Ritsch 
Institut fur Theoretische Physik, Universitdt Innsbruck, Technikerstrasse 25, A-6020 Innsbruck, Austria 

(Dated: 29 October, 2007) 

We study the effect of defects in the Mott insulator phase of ultracold atoms in an optical lattice 
on the dynamics of resonant excitations. Defects, which can either be empty sites in a Mott insulator 
state with one atom per site or a singly occupied site for a filling factor two, change the dynamics 
of Frenkel excitons and cavity polaritons. While the vacancies in first case behave like hard sphere 
scatters for excitons, singly occupied sites in the latter case can lead to attractive or repulsive 
scattering potentials. We suggest cavity polaritons as observation tool of such defects, and show 
how the scattering can be controlled in changing the exciton-photon detuning. In the case of 
asymmetric optical lattice sites we present how the scattering effective potential can be detuned by 
the cavity photon polarization direction, with the possibility of a crossover from a repulsive into an 
attractive potential. 

PACS numbers: 42.50.-p, 71.36.+C, 71.35.Lk 



I. INTRODUCTION 

In one of the most significant achievements of cold 
atom physics in recent years, a degenerate gas of ultra- 
cold atoms loaded into an optical lattice formed far off 
resonance lasers was demonstrated to undergo a quantum 
phase transition from a superfluid into the Mott insulator 
phase [l[ . In the Mott insulator a perfectly regular lattice 
at half the laser wavelength with a fixed number of atoms 
per site is achieved. This quantum phase transition from 
the superfluid into the Mott insulator was predicted and 
described by the Bose-Hubbard model As the Mott 

insulator phase can be considered as an artificial crystal 
with parameters well controllable in time and space, this 
generates a close connection between cold atom optics 
and solid state physics, which has generated a plethora 
of subsequent work generalizing the model and identify- 
ing more and more solid state phenomena to be studied 
in such configurations 0]. Furthermore, strong resonant 
light-matter coupling for a Bose-Einstein Condensate in 
an optical cavity has been achieved recently || . 

In previous work Q we exhibited the similarities of 
such artificial crystals and molecular or noble atom crys- 
tals, where one gets the so called Frenkel-excitons. In 
particular for the special case of identical ground and ex- 
cited state optical lattice potentials the Mott insulator 
phase allows the formation of such excitons as collective 
electronic excitations in the whole optical lattice @, HI • If 
the optical lattice is placed between cavity mirrors with a 
single cavity mode close to resonance with the excitons, 
coherent superpositions of excitons and cavity photons 
can form. In the strong coupling regime the correspond- 
ing system eigenstates then are cavity polaritons, which 
have very similar behavior in solid state and cold atom 
systems 0, HJ. Going beyond the analogous solid state 
case, very interesting new physics appears in the case of 
lattices with two atoms per site, where the on-site res- 
onant excitations form symmetric (bright) and antisym- 
metric (dark) superposition states of individual excita- 
tions [lfj. While the symmetric states radiate and via 



dipole-dipolc coupling will form excitons and thus cav- 
ity polaritons, the antisymmetric states have no dipolc 
moment, so that their excitations are long lived and stay 
localized at one site. Lattice asymmetries can lead to 
coupling between bright and dark states and produce 
strongly polarization dependent spectra. 

Optical lattices are spatially and externally described 
by the laser field and thus almost perfect over long dis- 
tances. Nevertheless due to imperfections in the dynam- 
ical formation of the Mott insulator, the appearance of 
some defects in the Mott insulator phase is unavoidable 
[l[ and there are missing or extra atoms at some sites. 
Such defects might create decisive errors when using such 
lattices for quantum information processing or for observ- 
ing even more exotic quantum phases in such systems and 
some suggestion to repair such defects were already pre- 
sented 111]. On the other hand they are hard to detect 
by direct observations of the atomic distribution. 

In this paper we investigate the effect of such defects 
on the dynamical properties of excitons and cavity polari- 
tons. We concentrate in the case of a low defect number, 
where the exciton and cavity polariton picture still holds 
and these quasi-particles are only scattered by single de- 
fects. Namely, the distance between each two defects is 
enough large for the formation of coherent excitons and 
cavity polaritons, which propagate as free quasi-particles 
between each two scattering processes. In principle de- 
fects can move by hopping among the lattice sites, but 
the corresponding time scale is so long that they can be 
considered frozen. Introducing more and more defects 
with faster hopping would then correspond to the tran- 
sition from the Mott to a superfluid atomic phase in the 
lattice, which will not be considered here. 

Here we treat defects related to a missing atom: in the 
case of one atom per site this corresponds to vacant sites 
and while for two atoms per site a defect is a singly oc- 
cupied site. Interestingly in calculating the exciton and 
cavity polariton elastic scattering amplitude of such de- 
fects, the two case behave even qualitatively very differ- 
ent. As a consequence of this defect scattering we show 
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how cavity polaritons can be used as observation tool to 
detect defects in the Mott insulator phase. 

The paper is organized as follows. In section 2 we in- 
vestigate the exciton scattering off a vacancy in an optical 
lattice in the Mott insulator with filling factor one, which 
is generalized to the scattering of cavity polaritons in sec- 
tion 3. The scattering of excitons and cavity polaritons 
off a defect in the Mott insulator case with two atoms 
per site is treated in section 4, where also we treat the 
special case of an optical lattice with asymmetric sites. 



FIG. 1: An impurity in ID optical lattice. 



II. EXCITON SCATTERING OF A VACANCY 
IN AN OPTICAL LATTICE 

Let us consider first a 2D optical lattice at the magic 
wavelength filled with one 2-level atom per site in the 
Mott-insulator phase. Including dipole-dipole interac- 
tions resonant excitations of the atoms can be repre- 
sented as quasi-particles called Frenkel excitons @. A 
single atom missing at the origin (r^ = 0) then creates 
an impurity in the artificial lattice of ultracold atoms, 
which for sufficient lattice depth stays localized and will 
not hop among the lattice sites as shown in figure (1). 
We will now consider scattering of these excitons at such 
a missing atom (hole). As the scattering time is much 
shorter than the hopping time, the impurity assumed to 
be localized during the scattering process. 

Adding the impurity to the ideal atom crystal just 
needs a small change in the system Hamiltonian H = 
Hq + V derived in Ref. |6j , whose components now read: 

H a = toex(k) BtB k , V = -hu A BlBo. (1) 



Hq represents free excitons due to the resonant electronic 
excitation transfer among lattice sites induced by dipole- 
dipole interactions including the lattice symmetry. 
and -Bk are the creation and annihilation operators of 
an exciton with in-plane wave vector k, respectively, and 
uj ex (k) is the exciton dispersion. V is the impurity Hamil- 
tonian at the origin, where Tuva is the atomic transition 
energy. 

For an exciton the impurity thus appears as a potential 
well located at the origin, with depth Hloa and radius a. 
Here we will not try to find the self-consistent eigenstates 
of H but to consider only the scattering problem. We can 
neglect trapping of an exciton in the impurity potential 
as a vacant site cannot absorb the trapping energy. Hence 
we consider only the scattering process of an exciton off 
the impurity. 

In the following we calculate the scattering amplitude 
of an exciton of an impurity. The initial exciton has wave 
vector k very far from the impurity, is scattered elasti- 
cally and found with wave vector k' also very far, and 
for elastic scattering we have |k'| = |k|. As the impu- 
rity is very deep with depth Hloa, in order to calculate 
the scattering amplitude we need to go beyond the Born 



approximation [12j even though the perturbation is con- 
fined to a radius a. 

The eigenstates of the free exciton Hamiltonian are 
Ho\4>k) — E{k)\(j)\^), and the full Hamiltonian eigenstates 
obey H\ijj) = E\tp). The excitation eigenstates in the 
quasi- momentum space |</>k) are related to the lattice 
space eigenstates by = -^J2i elkri where N 
is the number of lattice sites. The eigenstates form a 
complete basis, with the closure relation J^i \<f>i){<l>i\ = 1- 
A general state \ip) can be expanded by = £\ tpi \4n), 
with the expansion amplitude tpi = {4>i\^). 

In the scattering problem the incoming exciton is pre- 
pared in a delocalized unperturbed exciton state |</>k) 
and the scattered exciton is observed also very far from 
the impurity. Hence, the whole system eigenstates 
needs to obey the boundary condition \ip) — > |0k) as 
V — > 0. The required solution is given by the Schwinger- 
Lippmann equation fl2| 



where 



Gq = lim 

7)^o+ \ E - Ho + ir) 



(2) 



(3) 



The sign of +ir] is chosen in such a way to ensure that the 
scattered exciton propagates away from the impurity. In 
multiplying the equation from the left by the bra vector 
(4>i | , and in inserting the above identity operator between 
Go and V, we get 



E 

3 



{^\Go\^){^\V\^). (4) 



The calculations using the above definitions yields: 
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E, 



N 



N 



y lim 



e 



E-E(k') 



IT) 



V'o, (5) 



where E a = Hl>a- This represents the scattered wave 
amplitude at site i very far from the impurity. The first 
term is for the non-scattered wave, and the second is for 
the scattered wave. For the impurity site amplitude, at 
r = 0, we get 




N E r/™ + E - E(k') + ir] 



(6) 
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As the scattering is elastic the energy E can be replaced 
by the incident exciton energy and we have E — E{k). 
We get 



E a I 



(7) 



'N I l-E a I st 
where we defined the oscillating and static summations 



L, = lim 

k' 

ist = -T?y^ um 



ik'- 



+ E{k) - E{k') + if] 1 



1 



E{k) - E(k') 



17] 



(8) 



Our main aim now is to calculate the two sums in the case 
of small wave vector scattering excitons, that is ka <C 1. 
For this we can use the parabolic approximation for the 



dispersion of the excitons, u> ex {k) = u ex (0)- 



2m. 



where 



m ex is the exciton effective mass. For isotropic atoms in 
a square lattice of cubic symmetry, and in including only 
nearest neighbor interactions with coupling parameter J, 
we get m ex = —h/(2Ja 2 ), and ui ex (0) = oja — 2 J. For a 
large enough 2D optical lattice, the sum over k' can be 
approximated by a 2D integral given by 



N 4^ V2tt 



d 2 k'. 



(9) 



We get the two integrals 



Ios 



2m 



f-V 

h 2 \2itJ 



h 2 

2m, 



d 2 k' lim 
d 2 k' lim 



j)-»o + fc 2 - k' 2 + ir] ' 



as shown in the appendix. The scattered exciton, follow- 
ing the definition in [l3j . is given by 



1>i 



1 



N 



f(k) 



2fcr 



Akr 



(ii) 



where r = |r,| is the distance of the i site from the origin. 
We defined the scattering amplitude by 



/(*) = 



?rE a 
2A, , 



and the effective exciton band width is defined by 



h 2 ir 2 
2m ex a 2 



(12) 



(13) 



The dipole-dipole interaction energy between differ- 
ent sites in optical lattices is small, and hence the ex- 
citons have a small band width. In our case we have the 
limit of E a 3> A. ex so that the scattering amplitude is 



f(k) 



, which exactly reproduces the result for 



the scattering off a hard disk of radius a [14[, with the 
scattering cross section defined by a(k) = 27r|/(fc)| 2 . 

We thus conclude that an impurity generated by a 
missing atom in an optical lattice acts effectively like 
a hard disk of radius a. As the scattering is elastic, 
for an incident exciton with a fixed and small wave vec- 
tor k we get a scattered exciton with a wave vector k', 
which equals in magnitude to the incident one, namely 
|k| = |k'| = fc. Hence if a big number of incident excitons 
with identical wave vector are scattered off the impurity, 
we get a ring of radius fc of scattered excitons. 

The scattering amplitude is plotted in figure (2) as 
a function of wave vectors, fc, for lattice constant a — 
2000 A. The singularity at fc = is the 2D signature. 
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FIG. 2: The scattering amplitude vs. wave vector k, for zero 
detuning. 



III. CAVITY POLARITON SCATTERING OF A 
VACANCY IN AN OPTICAL LATTICE 

Let us now add cavity mirrors to our lattice as de- 
scribed in @, where in the strong coupling regime, the 
system eigenstates are cavity polaritons. As polariton is 
a coherent superposition of an exciton and a photon it 
can scatter of the impurity due its excitonic part. The 
corresponding Hamiltonians then read: 

H = J2 «M*0 4r^kr, V = -E a BlBo, (14) 

kr 

where Q±(k) are the upper and lower polariton branch 
dispersions, given by: fl±(k) = " c ""( fc )+" e:c ( fc ) ±Afc. Here 
Afc = y% + l5fc| 2 , and the exciton-photon detuning is 
defined by 5k = " c °"( fc )~""( fc ) , The parameter is for 
the exciton-photon coupling, which is taken to be of the 
electric dipole interaction. The cavity-photon dispersion 



,{k) 



where L is the 



is given by uj cavy n,, — ^= y ,„ , KL/ 
distance between the cavity mirrors, and e is the cavity 
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medium dielectric constant, which here taken to be a vac- 
uum with e = l. The polariton operators are defined by 
A^± = Xj^ Bk + Y^ 1 ak, with at the cavity photon oper- 
ator. X^ and are the exciton and photon amplitude, 

respectively, which are given by X^ = ±\J A £^*^ , and 

~ 2A fc (A fc =F^fc) ( more details are found in Ref.[6|). 
As next step we now calculate the polariton scattering 
off such impurity. Again the scattering is elastic so that 
an incident polariton in branch r with wave vector k will 
scatter into a polariton with wave vector k' in the same 
branch. 

The free eigenstates obey Hq\4>^) — E r (k)\<f^.), 
and the whole system eigenstates obey H\ip) — E\ip), 
where E r (k) = Ml r (k). The cavity photon eigenstate 
is defined by a\\vac) = |0k™), and the exciton eigen- 
state in quasi- momentum space is defined by B^\vac) — 



and the sums converted to the integrals 



The r polariton branch eigenstate is defined by 
= |<^?° ), and in terms of exciton and photon 



6-) 



states we get = X^ |< 

real lattice space, we get \4>f° l 



Y, 



± 



6£ au ). In terms of 



A'f 



k ± 7 y/N Ei 

Y± \<j)^ v ). A general state can be expanded as 

As above the scattered polariton state is given by the 
Schwinger-Lippmann equation — \cf%° r ) + GoV\i/j), 
where we used the free Green function operator in Eq.([3]). 
We multiply the equation from the left by the lattice 
exciton eigenstate {4>l x \, and insert the previous identity 
operator between Go and V, to get 

(trw = (^ric') +E^riG ic)(Ciw- (15) 



Using the above definitions, we get the scattered polari- 
ton excitonic part amplitude by 



^ = ^L e ik., i _Eay lim 



\X s y\ 2 e 



2 „ik' ri 



N 



N r^o + E r (k) ~ E s {k') + ir\ 



where due to energy conservation we replaced the energy 
E by the initial energy E r (k), and \Xf., | 2 is the exciton 
weight in the polariton. For the impurity site amplitude 
we have 



^ = J 1 + || y Um -, , ^ k 
y/N\ N frv-+o+ E r {k)-E, 



We thus obtain 




^ex _ X'k j ik-n , E a I os 

y/Nl 1 - E a I st 



(17) 



(18) 



where we defined oscillating and static sums. For a large 
optical lattice size the fc-space can be assumed continuous 



Ins 



1st 



\X°,\ 2 e*'- 



-Y(-) / d 2 k' lim 
^\2tt/ J 7,^0+ E r {k) - E s (k') + ir) 



d 2 k' lim 



v->o+ E r (k) - E s {k') + ir]' 
(19) 



We calculate first the oscillating integral in case of scat- 
tering of small wave vector polaritons, that is in the limit 
ka <C 1 in the lower branch. In this limit the polari- 
ton dispersion can considered approximately parabolic 
with a polariton effective mass of m p . This is of the or- 
der of the cavity photon effective mass m p rs (Tuk)/{cL). 
Hence, the lower polariton branch dispersion is taken to 
be E-(k) = S_(0) + 1^-. The scattered states will also 
have small wave vectors due to the energy conservation. 
As the upper branch has higher energies its contributions 
are negligibly small here. Anyway the upper branch for 
larger wave vectors is mainly photonic with small exci- 
tonic part and thus only weakly contributes to the scat- 
tering of impurities. We consider the case around zero 
detuning between the excitons and photons. Therefore 
the excitonic weight |AjT| 2 will change from half around 
zero wave vector up to one for large wave vectors, where 
the lower branch became excitonic. The weight |A^| 2 is 
a smooth function of k, and in the present case can be 
taken out of the integral, and to be fixed with the ini- 
tial polariton excitonic weight value. From here on we 
neglect the contribution of the upper polariton branch, 
both as initial and scattered states, and also as interme- 
diate scattering state. We will drop the branch index and 
summation, and all the parameters will be only for the 
lower branch. 

The integral now reads: 



I or — —X, 



2 ^ifTLip 



d 2 k' lim 



ik'-: 



7j^o + k 2 — k' 2 + it] 



(20) 



This is similar to the case of excitons of the previous 
section except for the factor | A£| 2 and involving the much 
smaller polariton effective mass in place of the exciton 
mass. Hence we have m p -C m ex and the integration 
(see the appendix) gives 



r -_**2 
° s ~ 2A„ 



2kr 



■ikr 



(21) 



where we defined the polariton effective band width A p = 

h IT 

2m p a 2 ' 

We turn now to the second static integral. Here no 
oscillating exponent appears, and hence all wave vectors 
along the lower branch contribute, so that the interme- 
diate scattering state can be anywhere along the lower 
branch. Again for energy reasons we neglect the con- 
tribution of the upper branch. The integral then can 
be simplified by using a model for the lower branch dis- 
persion in place of the real one [l5j . The lower branch is 



5 



divided into two parts. The first part between < k < ko 
is taken to be of a parabolic dispersion with a polariton 
effective mass, where E(k) = E(0) + (H 2 k 2 )/(2m p ); and 
the second part between ko < k < ir/a, where ir/a is the 
Brillouin boundary, is taken to be dispersion less with 
energy equal to the exciton energy at zero wave vector, 
where we have feo «C n/a. The intersection point between 
the two parts is fixed by h 2 k 2 /(2m p ) — E — E(0). The 
integration as shown in the appendix gives 



T -_**2 





7T 


Hi) 


— 1 — 




2 



3a? < 22 > 



where A& = Eq —E(k). The first term is the contribution 
of the parabolic part and the second originates from the 
flat part. The Ln(k/ko) is large in the limit ka <C 1, but 
for a laterally confined 2D optical lattice, and for the 
smallest wave vector k m which is different from zero, due 

to the fact A p ^> E a , the term 2 £ m In 

smaller than one and can be neglected. 
7r/(4A fe ). 

The scattered polariton excitonic part reads 



(jho) ^ S mucn 
Hence I st w 




(23) 



where the scattering amplitude is given by 



/(*) 



tyE„ 
2A, 



1 - 



irE a 
4A fc 



(24) 



k ~ polaritons become more photonic and they scat- 
tered off the impurity much less. As the atomic energy E a 
is much larger than the coupling energy hg, the scatter- 
ing amplitude only shows a negligibly small asymmetry 
for negative and positive detunings, but we expect a clear 
resonance of the scattering amplitude around resonance 
8q. Considering a stream of incident polaritons with the 
same wave vector k = fck we then should see a ring of 
radius k of scattered polaritons. 
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FIG. 3: The scattering amplitude vs. exciton-photon detun- 
ing, for k w polaritons. 



As E a > A fc , we have (nE a ) j (4A&) > 1, and the scat- 
tering amplitude is negative, that is f(k) < 0. Again 
the above impurity effective potential is repulsive for the 
polaritons. 

In the case of zero detuning, that is 5o = 0, we have 
for small wave vectors X 2 = 1/2, and Afc is of the order 
of the exciton-photon coupling, \gk\, where E a 3> Afc, 
and we get f(k) » — Afc/A p . This scattering amplitude 
equals that of the scattering from an effective potential 
of a square barrier potential of height Afc and width a. 

For the case of negative detuning, in the limit of E a ^> 
<5o, where Afc is of the order of the exciton-photon detun- 
ing, the scattering amplitude is f(k) — (2X 2 A k )/A p . 
This scattering amplitude equals that of the scattering 
from an effective potential of a square barrier potential 
of height 2X 2 A k and width a. Hence we find that the 
scattering amplitude can be controlled by changing the 
exciton-photon detuning. 

Lets insert some numbers: for an atomic energy 
E a = 2 eV, exciton-photon coupling H\g\ — 0.0001 eV, 
and lattice constant a = 2000 A, the scattering am- 
plitude is plotted in figure (3) as a function of the 
exciton-photon detuning, for approximately normal in- 
cident waves, where we used k = 10~ 6 A. It is clear 
that maximum scattering is obtained for zero detuning. 
As the detuning increases, both positive or negative, the 



IV. EXCITON AND POLARITON 
SCATTERING OF A VACANCY IN TWO-ATOMS 
PER SITE OPTICAL LATTICES 

We will now go one step further and study the scatter- 
ing of excitons and polaritons in a doubly occupied opti- 
cal lattice where one atom is missing in one site (called 
the origin site). The corresponding exciton polariton pic- 
ture for an ideal lattice is discussed in Ref.[l(J. Here 
resonant electronic excitation transfer occurs among the 
two on-site atoms with transfer parameter Jo and also 
among the nearest neighbor sites with transfer param- 
eter J\. The on-site interaction can be diagonalized to 
form symmetric and antisymmetric localized excitation 
with frequencies uj s = u>a + Jq, and uj a — u>a ~ Jo, re- 
spectively. Only the symmetric states have a significant 
dipole moment and thus form excitons through resonant 
coupling. The symmetric exciton dispersion in the limit 
ka <C 1 then reads to s (k) = uja + Jo + 8Ji + 2i ^ fc ■ , with 

the effective mass m e // = —h/ (iJia 2 ) and up to a rescal- 
ing is analogous to the previous case. Naturally within 
a cavity only the symmetric excitons are coupled to the 
photons to form cavity polaritons (for details see [HI). 
At first sight there seems to be nothing new in this case. 
However, we will show that the defect, plotted schemati- 
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cally in figure (4), which in this case constitutes a single 
atom as compared to two, here acts substantially differ- 
ent from the simple no atom hole in the case of a single 
occupied lattice. Physically this can be understood from 
the fact that a single atom can store almost the same 
amount of energy as a pair if only one photon excitations 
are considered. 

Let us now calculate the scattering amplitude of exci- 
tons and polaritons of this impurity. The scattering is 
for the symmetric excitons, and for the polariton exci- 
tonic part. The antisymmetric states are not involved in 
the scattering process, as they are on-site localized and 
posses no dipole moment. 



Antisymmetric state 
X^Symmetric state 
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FIG. 4: An impurity in ID optical lattice. 

Again we concentrate on the scattering of long wave 
length excitons with a parabolic dispersion described by 
an effective mass. The scattering amplitude is found to 
be 



/(*) 



7T Jo 

2A, , 



1 



2A ca 



[In I 



(25) 



For negative on-site energy transfer, Jo < 0, the scat- 
tering is of an effective potential of a square barrier of 
height | Jo I and width a. One need to compare between 
the on-site transfer parameter \Jq\ and the exciton band 
width A ex . In the limit of |J | ^> A ea . the scattering is 
identical to that from an effective potential of a hard disk 
with radius a (see figure (2)), that is f(k) « in ^ ka ^ ■ 

For polaritons the scattering of the impurity is only for 
their excitonic part. We consider the scattering of long 
wave length lower branch polaritons with a parabolic dis- 
persion and an effective polariton mass m p . The scatter- 
ing amplitude is given by 



/(*) 



IT Jo 

2A V 



1 _ [ Eh. 

1 ' 4A fc 



(26) 



Here Afc = Eq — E(k), where Eq = Hda + hJo, an d where 
Jo > Ji- 

In the limit of |J | > A k we get f(k) w -X% 

where, for Jo < 0, the scattering acts like a square barrier 
of height 2X?Afe and width a. Similarly in the opposite 



limit of | J 1 < A k we get f(k) « X\ (f^), and the 
scattering, for Jo < mimics a square barrier of height 
7r| Jq\X^/2 and width a. As shown in the previous chap- 
ter, the polariton scattering of an impurity can be modu- 
lated by controlling the exciton-photon detuning. As we 
change 8k, then X\ and A& are strongly changed, and as 
a result the scattering amplitude f(k) and the effective 
potential strength are changed. 

Introducing the numerical examples used in the pre- 
vious section with transfer energy HJq = —0.001 eV the 
scattering amplitude is plotted in figure (5) as a func- 
tion of the exciton-photon detuning, for incident waves, 
namely k s» 0. It is clear that maximum scattering is ob- 
tained for zero detuning. As here the transfer energy HJq 
is smaller than, and close to, the coupling energy hg, the 
scattering amplitude asymmetry for negative and posi- 
tive detuning is much more pronounced relative the case 
of the previous section. For k = 0, at zero detuning, the 
lower branch polariton is half exciton and half photon. 
For negative detuning, at k = 0, the lower branch polari- 
ton is more photonic than excitonic; while for positive 
detuning the opposite the lower branch polariton is more 
excitonic than photonic. Therefore, the positive detuning 
scattering amplitude is larger than the negative one. 
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FIG. 5: The scattering amplitude vs. exciton-photon detun- 
ing, for k m polaritons. 



A. An impurity in an asymmetric optical lattice 

One of the most interesting phenomena for excitons in 
doubly occupied lattice s ap pears in the case of asymmet- 
ric optical lattice sites [1J| . We get a strong polarization 
dependence of the exciton and polariton energies. Here 
we discuss the consequence of this for their scattering 
amplitudes off impurities. 

If one of the orthogonal pairs of counter propagating 
lasers forming the lattice potentials has a different inten- 
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sity, the potential is elongated in one direction, e.g. x. 
Hence the two atoms at one site will have a larger average 
distance R in the x direction as shown in figure (6). The 
atomic transition dipole induced by the cavity photon 
is assumed to be in the (x,y)-plane fl — /i(cos 9, sin 6), 
where 9 is the angle between the dipole /2 and the x 
axis. The resonance dipole-dipole interaction between 
the two on-site atoms is HJq{6) — HJ (1 — 3 cos 2 9), 

-2 

where hJ = R3 . We have HJ (9 = 0) = —2ft J, and 
HJ (9 = 90) = HJ, with hJ {9 w 54.74) = 0. The detun- 
ing energy between symmetric and antisymmetric states 
thus will change sign at some angle, where the polariza- 
tion even vanishes. 
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FIG. 6: An asymmetric optical site. 

The long wave length polariton scattering amplitude 
as a function of the relative transition dipole direction 
now is 

fk{ 0) = '} 2A \) . (27) 
\4A k (e)J 

We plot this scattering amplitude as a function of 9 for 
different values of J, and for k « 0, again using the 
numerical values of section 3. We also have E c (0) = 
Ea+J{0), and h,] 1 = 1(T 7 eV. In figure (7) the plot is_for 
hJ — 10 -4 eV, where h\g\ = hJ. Figure (8) is for hJ = 
5 x 10 -4 eV, where h\g\ < HJ. Note that the scattering 
amplitude changes sign at 9 m 54.74 and we can turn of 
scattering by using this proper angle. In figure (9) we 
plot this for HJ = 10~ 3 eV. Here two resonances appear. 
Figure (10) is for HJ = 5 x 10 -3 eV. Now the resonances 
tend to the negative-positive crossover angle. 

The negative scattering amplitude corresponds to a re- 
pulsive effective potential, and polaritons scattered away 
of the impurity. While the positive scattering amplitude 
corresponds to an attractive potential. Here even the for- 
mation of bound states could be expected, allowing for 
localized polaritons. The bound state signature appears 
in the pole of the scattering amplitude of Eq. (f2"T)) but by 
including the Ln term of Eq. lffi?]) in the denominator. We 
get a shallow bound state that falls inside the polariton 
line width, and which presents as a resonance scattering 
state inside a continuum. 



V. SUMMARY 

We demonstrated that defects in an optical lattice in 
the Mott insulator phase will change the dynamics of 
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FIG. 7: The scattering amplitude vs. 9, for k ~ polaritons, 
with hJ = 10" 4 eV. 
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FIG. 8: The scattering amplitude vs. 9, for k ~ polaritons, 
with fij=5x 1(T 4 eV. 



resonant excitons and cavity polaritons. For a very low 
density of defects, the exciton and polariton pictures still 
hold and the defect simply acts as a scatterer for such 
quasi-particles. We calculated the scattering amplitude 
for excitons and polaritons showing that the defects can 
be approximated by effective potentials. In the case of 
one atom per site the vacancy scattering effective poten- 
tial for excitons is a hard disk of a radius equal to the 
lattice constant a, which is quite different for polaritons, 
where only the excitonic part contributes significantly 
to the scattering. The scattering effective potential for 
long wave length polaritons just is a square barrier of 
height 2|Xfc| 2 Afc and width a, where \Xk\ 2 is the exciton 
weight in the cavity polariton, and which appears here 
as the scattering is only for the polariton excitonic part. 
For the case of zero exciton-photon detuning, and in the 
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FIG. 9: The scattering amplitude vs. 9, for k ~ polaritons, 
with hJ = 1(T 3 eV. 
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FIG. 10: The scattering amplitude vs. 8, for fe « polaritons, 
with SJ = 5x 10" 3 eV. 



presented for the case of asymmetric optical lattice sites. 
Also we showed that at a fixed angle Jo can change sign 
from positive to negative, and then also the scattering 
amplitude changes sign. 

We conclude that cavity polaritons can be used as a 
useful tool to observe defects in an optical lattice. As 
the polariton is part exciton and part photon, where the 
photon part can be controlled externally, the above re- 
sults suggest that one can recognize the kind of defect in 
each case. For example in transmission or reflection ex- 
periments, for an incident held with a fixed wave vector 
k, in the transmitted or reflected signals we get a scatter- 
ing ring of radius k with the appropriate intensity which 
obtained from the scattering cross section. 
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APPENDIX A: THE CALCULATION OF EQ. (flO|) 
INTEGRALS 



Here we give the calculation details of the two integrals 
of Eqs. lfTOl) . 

The oscillating integral is written explicitly as 



2m e; 



(£)" 



+ 00 



dk' 



d6 lim 



j^/gifc'r sin t 



rj— >0+ k 2 — k' 2 + IT] ' 

(Al) 

where r = | | is the distance of the i site from the ori- 
gin, and k' — |k'|. The integral over k' is from zero to 
infinity, which is possible here due to the oscillating ex- 
ponents, where the contributions of large wave vectors 
are negligible. The integral over 6 is calculated by 



Jo(k'r) = 



1 

2^ 



d6 



Ah r 



(A2) 



limit ka <C 1, we have \Xk\ 2 = 1/2, and A& is of the 
order of the Rabi splitting. While for large negative de- 
tuning we have \Xk\ 2 <C 1, and A& is of the order of the 
exciton-photon detuning. As \Xk\ 2 and A& are a function 
of the detuning, hence the scattering amplitude can be 
controlled by changing the detuning, with a scattering 
resonance at zero detuning. 

In the case of a two atoms per site lattice with a defect 
represented by a single atom site, we found two distinct 
cases. If the on site dipole-dipole coupling Jo is larger 
than the Rabi splitting, we got similar results as in the 
case of one atom per site. But if Jo is smaller than the 
Rabi splitting, hence the polariton scattering effective 
potential is identical to the scattering of a square barrier 
of height Jo and width a. The parameter Jo can be eas- 
ily controlled in changing the photon polarization, as we 



where Jo(x) is the zero order Bessel function of the first 
kind. We have now 



In 



2m e; 



dk' lim 



k'Jo(k'r) 



We use the result 
dk' lim 



fc , J (fc , r) 
r\— >o + k' 2 — k 2 — irj 



»;^o + k' 2 — k 2 — irj 

(A3) 



l |^ 1) (fcr), (A4) 



where Hfi'(x) is the zero order of the Hankel function 

of the first kind, which is defined by H^\x) = Jq(x) + 
iNo(x), where Nq(x) is the zero order Bessel function of 
the second kind. The scattered wave is observed very far 
from the impurity, that is at r — > +oo. For large x we 
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have the asymptotic expansions 
~2~ 



J (X) r 

which yield 
The result is 



.(x - |) , N (x) 



J(x- ir/4) 



2A e 



7.7T 

2fcr 



Akr 



dn (z - J) , 
(A5) 

(A6) 



(A7) 



where we defined the effective exciton band width in 
Eq.®. 



For the static integral the angle integral gives 2-7T, we 
have 



1st = 



2m eL 



7r I a 



dk' lim 



v^o + k' 2 



(A8) 

Here, due to the fact that no oscillations exist, it is not 
possible to extend the integral to infinity. We used a 
cut-off at the Brillouin zone boundary. Such cut-off is a 
physical one as the maximum exciton wave vector is for 
wave length equals the lattice constant. The integration, 
in the limit of ka <C 1, yields 



2 A» 





7T 











(A9) 
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